Measurement of gravitational spin-orbit coupling in a binary pulsar system 
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In relativistic gravity, a spinning pulsar will precess as it orbits a compact companion star. We 
have measured the effect of such precession on the average shape and polarization of the radiation 
from PSR B1534+12. We have also detected, with limited precision, special-relativistic aberration of 
the revolving pulsar beam due to orbital motion. Our observations fix the system geometry, including 
the misalignment between the spin and orbital angular momenta, and yield a measurement of the 
precession timescale consistent with the predictions of General Relativity. 


PACS numbers: 04.80.Cc,97.60.Gb,97.80.Fk,95.85.Bh 

Parallel transport of the angular momentum of a gyro- 
scope moving in curved spacetime leads to geodetic pre- 
cession [1]. In the solar system, the only observed exam- 
ple is precession of the Earth-Moon system as it orbits 
the Sun [2]. The recently launched Gravity Probe B ex- 
periment plans to measure the geodetic precession of a 
gyroscope in Earth orbit as well as its gravitomagnetic 
Lense-Thirring precession [3]. 

Binary pulsar systems are important laboratories for 
gravitational physics [4], in part because the strong self- 
gravity of neutron stars ( GM/Rc 2 ~ 0.2) raises the pos- 
sibility of deviations as large as order unity in some al- 
ternate gravity theories — even theories that agree with 
general relativity (GR) in weak-field tests [5]. Imme- 
diately after the discovery of the first binary pulsar, 
PSR B1913+16 [6], it was realized that geodetic preces- 
sion could lead to variations in the path of the observer’s 
line of sight across the pulsar’s magnetic pole, and hence 
changes in the radiation pattern at Earth [7]. Recently, 
the anticipated' pulse-profile variations have been ob- 
served, in good qualitative agreement with predictions 
[8-10], but uncertainties in the intrinsic beam shape pre- 
vent a quantitative measurement of the precession rate. 

Here we describe observations of the double-neutron- 
star system PSR B1534+12, in which precession is chang- 
ing the observed pulsar profile by about 1% per year 
[11, 12]. We outline a new, general technique for com- 
bining measurements of long-term quasisecular varia- 
tions caused by geodetic precession with those of pe- 
riodic variations induced by special relativistic aberra- 
tion modulated by orbital motion. Together, the obser- 
vations allow a quantitative estimate of the precession 
rate, independent of the unknown pulsar beam shape. 
Although the precision is still low, we show that this 
model-independent precession rate is consistent with the 
predicted rate in GR. The polarization properties of the 
pulsar signal are also changing: with the additional as- 
sumption that the magnetic field structure is dipolar, we 
can determine the angles between the spin and orbital 



FIG. 1: The spin and orbited geometry we derive for 

PSR B1534+12, after [4]. The pulsar spin axis is S and 
the orbital angular momentum direction is k, with misalign- 
ment angle 6 (not shown). The angle between S and the 
line of sight is the supplemental angle A = 7r — £ is 
shown. The angle between S and the magnetic pole is a, 
and (3 = £ — a is the minimum impact angle of the magnetic 
pole on the line of sight. The projection of S on the plane 
of the sky provides angle tj, measured counterclockwise from 
the ascending node; —tj is shown. Spherical geometry gives 
cos 5 = — sin i sin A sin t] + cos A cos i. 

angular momenta and the line of sight to the pulsar. 

The spin and orbital geometry is shown in Fig. 1. The 
pulsar spin axis S will precess around the total angular 
momentum, which is very well approximated by the or- 
bital angular momentum direction, precession will cause 
potentially detectable periodic variations of the projec- 
tion of S on the plane of the sky, and on the inclination of 
S with respect to the observer. In GR, the time averaged 
precession rate of the pulsar can ,be written [7]: 

n spin _ 1 ( GMq \ 2/3 / P h \ ~ 5/3 m 2 (4 mi + Zm 2 ) 
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where G is Newton’s constant, M© is the mass of the Sun, 


c is the speed of light, mi and m 2 are the pulsar and com- 
panion masses, respectively, Pb is the orbital period, and 
e is the eccentricity. A description in generalized theories 
of gravity is given in [4]. For PSR B1534+12, using the 
stellar masses determined through high-precision timing 
[13], the precession rate predicted by GR is 0.51° /yr. 


OBSERVATIONS AND ANALYSIS 


Observations were made with the 300-m Arecibo radio 
telescope, using the “Mark IV” data acquisition system 
[14] at an observing frequency of 430 MHz. The signal 
was processed using coherent dedispersion, providing full 
polarization information as well as a pulse shape unaf- 
fected by dispersive smearing in the interstellar medium. 
Data acquisition details have been described elsewhere 
[13]- 

The data span the interval from mid-1998 to mid-2003, 
incorporating some 400 hours of observing time. The 
pulsar was observed biweekly or monthly. The observed 
signal strength varies widely because of interstellar scin- 
tillation; only epochs with high signal-to-noise ratio were 
used here. Campaigns of roughly 12 contiguous observ- 
ing days were also conducted every summer except 2002. 

The cumulative pulse profile from the representative 
2001 June epoch is shown in Fig. 2. Here we also show 
a fit of the linear polarization to the standard “rotating 
vector model” (RVM) [15], in which the position angle 
of linear polarization ip is assumed to be parallel to the 
plane of curvature of magnetic dipole field lines rotating 
with the star, giving 


tan [ip (0) - ip 0 \ = 


sin a sin ( 0 — 0o ) 


cos a sin ( - sin a cos C cos (0 - (po) ' 
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where 0 is the pulse phase and 0o and ipo are constants. 

The position angle sweep is observed over most of the 
pulsar period, and the model fit is generally good. There 
are strong deviations from the model near the pulse peak, 
as is often seen for “core” profile components [16]. We 
exclude -this region from our fits. The data are con- 
sistent with a roughly orthogonal rotator model, with 
a magnetic inclination angle a = 102.8 ± 0.5° and the 
line of sight passing between the magnetic pole and the 
stellar equator, within a few degrees of the magnetic 
pole. We also consider the time evolution of a and the 
impact parameter of the line of sight on the magnetic 
pole /?, using cumulative profiles from each campaign, 
one especially strong biweekly observation, and an ear- 
lier coherently-dedispersed profile from observations with 
the “Mark III” data acquisition system [17]. As expected, 
the data are consistent with no evolution of a. However, 
as shown in Fig. 2, is changing with time, at a rate 
dfi/dt = —0.21 ± 0.03°/yr, resulting in a larger impact 
parameter at later times. This change in (3 is direct ev- 



FIG. 2: Top panel: the position angle of linear polarization 
in 2001 June, measured clockwise on the plane of the sky (the 
convention in [4]), with best fit rotating vector model (RVM) 
overlaid. Only the position angle points indicated by large 
dots were used in the RVM fit; these were weighted by their 
uncertainties, with a small uncertainty added in quadrature 
to account for deviations from the RVM. Middle panel: total 
intensity (solid) and linear polarization (dashed) profiles in 
2001 June. This profile is very similar in shape to our “ref- 
erence” profile P 0 . Inset: evolution of impact angle p with 
time. Bottom panel: “Difference” profile Pi , representing es- 
sentially the time-derivative of the observed profile. 


idence of geodetic precession, and can be related to the 
system geometry and precession rate in a simple fashion 
[4]: dp/dt = D® pin cos tj sin i (see Fig. 1). 

The shape of the profile is also changing with time 
and orbital phase, allowing a completely independent 
probe of the precession. Secular changes in the profile 
were first noticed at 1400 MHz [12], but evolution of the 
430-MHz emission only became apparent with coherently 
dedispersed observations [11]. Shape variations are more 
difficult to connect directly to the precession rate than 
polarization changes. It can be done with an assumed 
model of the beam shape, as has been attempted for 
PSR B1913+16 [9, 10, 18], but in that case the true 
beam shape is still debated and the results are therefore 
less than satisfactory. 

Here we note that it is possible to make a model- 
independent precession estimate by also measuring the 
orbital modulation of the profile shape caused by 
aberration — a special-relativistic effect independent of 
strong-field gravity. Aberration shifts the observed an- 
gle between the line of sight and spin axis by an amount 
[4]: 


(JaC = [-cos TjSiu) + cos i sin p C (u)] , (3) 

sin i 
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where /?i = nxjy/{ 1 — e 2 ) is the characteristic velocity 
of the pulsar, with the orbital frequency n = 2i t/Pb, the 
projected semimajor axis x = ai sin i/c, and the eccen- 
tricity e all available from timing data, and where C (u) = 
cos[w± A e (u)]+e cos a; and S(u) = sin[a;±.<4 e (ti)]±esinu; 
are functions of the time-dependent angle of periastron 
passage u and the eccentric anomaly u through the true 

anomaly A e (u) = 2 arctan f ( 


\l-ej 


tan j 


Now let F(Q be any function defined on the observed 
pulsar signal that depends on the viewing angle (such as 
integrated intensity, component width, polarization frac- 
tion, etc.) . For small changes in the impact parameter, 
we Taylor-expand F(0 « F((o) + (F‘ where prime de- 
notes derivative with respect to £. The effects of aberra- 
tion and precession can then be written 


SaF 

dF 

dt 


F' [— cos T] S (u) + cos i sin n C («)] , (4) 
sin 7 

F'nf n sin % COS TJ. (5) 


The unknown beam shape enters only through F' , which 
can be eliminated by dividing these two equations. Mea- 
surements of both the orbital variation of F and its sec- 
ular drift thus allow tan 77 and fi* pm to be determined in 
a model independent way. 

For PSR B1534+12, we have measured the evolution 
of the totals intensity profile (Fig. 3). The strongest data 
scans from the annual campaigns were averaged into 12 
orbital phase bins, and analyzed together with the strong 
biweekly scans. Profiles with unusually low signal-to- 
noise. ratios or suspect calibrations were discarded. We 
used standard principal component (PC) analysis tech- 
niques [19] to derive orthogonal “reference” (Pg) and “dif- 
ference” (Pi) profiles that completely describe the pro- 
file evolution; Pi and a single-epoch profile very simi- 
lar to P 0 are shown in Fig. 2. The pulse profile P of 
each observation is well modeled as a linear combina- 
tion P = coPo + ci Pi. Because the overall amplitude 
varies, with scintillation, we chose as our observable quan- 
tity F the ratio Ci/co- The PC analysis provided es- 
timates of eg, ci, and their uncertainties, which we in- 
dependently checked through a frequency-domain cross- 
correlation technique in which a linear combination of 
the two profiles was fit in an iterative manner. Finally, 
we simulated the cross-correlation analysis to assess its 
sensitivity to systematic errors induced by imperfect cal- 
ibration or polarization cross-coupling, finding that such 
problems should be negligible in our dataset. 

The secular trend in F is evident in Fig. 3. This cor- 
responds to a decrease in the intensity of the core region 
of the profile relative to the lower level emission in the 
wings, consistent with early indications from 1400 MHz 
data, and as expected if precession is moving our line 
of sight away from the magnetic pole, as indicated by 
the polarization analysis above. The residuals after re- 
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FIG. 3: The shape parameter F = ci/co (see text) is shown 
as a function of date in the main panel and aberration phase 
(essentially the true anomaly corrected for the advance of pe- 
riastron) in the inset. The best-fit model is shown by the 
solid line in each panel, and in the orbital-phase plot, the GR 
prediction based on the RVM model is indicated by the dot- 
ted line. We have included a small error in quadrature with 
the measurement errors to account for smearing caused by 
averaging over a range of orbital phase. To account for sys- 
tematic errors, variable data quality, and unevfen time sam- 
pling, we have used a bootstrap analysis [19] to 7 estimate the 
uncertainties on model parameters. The resulting values and 
uncertainties are in good agreement with estimates obtained 
by scaling measurement errors to obtain a reduced-x 2 of 1. 


moving the best fit line axe shown as a function of orbital 
phase. A simultaneous linear fit of F as a function of date 
and of S{u) and C(u) gives the constraints flj Pm sin 2 i = 

0.42tg;i5 0 /yr and fl® pm sin i tan i cot 77 = OAitojj ° /yr 
( 68 % confidence), where we have used /?i = 0.67 x 10 -3 . 
A fit of F only as a function of date yields a x 2 value 
that is 15% higher than that for the full fit. 


DISCUSSION 

We have observed both long- and short-term varia- 
tions in the pulse shape of PSR B1534+12, as expected 
from geodetic precession and aberration. Assuming a 
dipolar field geometry and GR, the impact parameter 
change dfi/dt = f2® pm cos 77 sin i = —0.21 ± 0.03°/yr 
yields a measurement of the previously unknown an- 
gle 77 = ±115.0 ± 3.7°. Using the pulse timing value 
sin i = 0.975 [13] but making no assumptions about 
the validity of the RVM, our measured pulse profile 
variations yield, using equations 4 and 5, the consis- 
tent result 77 = ±103 ± 10°. Moreover, we may now 



solve for the precession rate, ftf 5111 = 0.44+°;jg o /yr (68% 

confidence) or = 0.44ig^ 3 °/yr (95% confidence). 
This value compares well with the GR-predicted rate of 
fl® pm = 0.51°/yr. 

The misalignment angle S between the spin and or- 
bital angular momenta can also be constrained. The 
angle A is known from the polarization studies. Only 
the absolute values of sin tj and cost are known, but our 
profile fit requires that cos i tan r] > 0. Therefore there 
are two possible geometries: i = 77.2° and r) = -115°, 
which gives <5 = 24.9 ± 3.8°, or i = 102.8° and r) = 115°, 
which gives S — 155.1 ± 3.8°. Both give identical results 
for precession in GR. As the angular momenta were al- 
most certainly aligned before the second supernova, the 
smaller misalignment value is favored on astrophysical 
grounds [20]. The preferred geometry then has i = 77.2°, 
r) = -115.0 ±3.7°, and <5 = 24.9 ±3.8° (Fig 1). The mis- 
alignment angle can be used to constrain mass loss and 
asymmetry in the second supernova. A full analysis will 
be published elsewhere, along with a study of the two- 
dimensional beam geometry of PSR B1534+12. 

Although the precision is, as yet, limited, this is the 
first beam-model-independent measurement of the pre- 
cession rate of a binary pulsar. It is also the first un- 
ambiguous reconstruction of the full three dimensional 
geometry of a binary pulsar system, as well as the 
first direct estimate of the misalignment angle. Future 
prospects for improvement include direct estimation of 
the angle 77 by combining scintillation studies [21] with 
polarimetry [4]. The derived geometry also allows us to 
predict the effects of aberration on the pulse timing [4]; 
this will in principle allow more precise timing tests of 
GR in future. 

We emphasize that the general technique of combining 
observations on the orbital and precessional timescales 
to make model independent precession rate estimates 
is potentially far more general than the particular ex- 
ample given here. An especially interesting prospect 
is the recently discovered highly relativistic system 
PSR J0737-3039 [22]; with /?i nearly twice as large as 
that of B1534±12, and a predicted precession timescale 
for the recycled pulsar of only 75 years, both effects will 
be quickly measured for this new system. 


* Electronic address: sta.irsQastro.ubc.ca; URL: http: 
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